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ABSTRACT 

We performed deep photometry of the central region of Galactic globular cluster M15 from archival 
Hubble Space Telescope data taken on the High Resolution Channel and Solar Blind Channel of the 
Advanced Camera for Surveys. Our data set consists of images in far-UV (FUVuo; F140LP), near- 
UV (NUV220; F220W), and blue (B 435 ; F435W) filters. The addition of an optical filter complements 
previous UV work on M15 by providing an additional constraint on the UV-bright stellar popula- 
tions. Using color-magnitude diagrams (CMDs) we identified several populations that arise from 
non-canonical evolution including candidate blue stragglers, extreme horizontal branch stars, blue 
hook stars (BHks), cataclysmic variables (CVs), and helium-core white dwarfs (He WDs). Due to 
preliminary identification of several He WD and BHk candidates, we add M15 as a cluster containing 
a He WD sequence and suggest it be included among clusters with a BHk population. 

We also investigated a subset of CV candidates that appear in the gap between the main sequence 
(MS) and WDs in FUVi 40 -NUV 2 2o but lie securely on the MS in NUV220-B435. These stars may 
represent a magnetic CV or detached WD-MS binary population. Additionally, we analyze our can- 
didate He WDs using model cooling sequences to estimate their masses and ages and investigate the 
plausibility of thin vs. thick hydrogen envelopes. Finally, we identify a class of UV-bright stars that 
lie between the horizontal branch and WD cooling sequences, a location not usually populated on 
cluster CMDs. We conclude these stars may be young, low-mass He WDs. 

Subject headings: globular clusters: individual (M15) - binaries: close - stars: Population II - white 
dwarfs - stellar dynamics - ultraviolet: stars 



1. INTRODUCTION 

The Galactic globular cluster (GC) MI 5 is the archety- 
pal core-collapsed GC with an extre mely small and dense 
core (as measured most recently by Ivan den Bosch et al.l 
2006) and thus has been an object of interest for many 
observers. Though convincing models for core collapse 
have existed for some time (e.g.. lCohr]|198C)IL many ques- 
tions remain and recent attention has turned to how core 
collapse might affect the stellar populations in the central 
regions of GCs. 

These high density central regions (MIS: po ~ 7 x 10 6 
M Q pc~ 3 : lvan den Bosch et al. 2006) provide an environ- 
ment that can be ideal for the production of many types 
of exotic binaries. In such a setting, the encounter rate 
is sufficiently high that one expects to find an increased 
number of close binary systems due to dynamical inter- 
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actions between stars that can har den primordial bina - 
ries and form new binary systems (jPoolev et al.l 12003). 
Because of the large amount of binding energy (some- 
times comparable to the binding energy of the cluster) 
contained in close binary systems, interactions of these 
systems can have a signi ficant effect on the dynamical 
evolution of the cluster (jHut et al.l 119921 ) and thus are 
of particular importance in understanding effects of core 
collapse. 

Close binary systems, especially those containing a de- 
generate remnant, are expected to be overabundant in 
GC cores due to both the preferential formation of close 
binaries in high density environments as well as the ef- 
fects of mass segregation. Many such systems are X-ray 
sources and therefore should be identified in X-ray sur- 
veys; M15 contains six X-ray sources with optical coun- 
terparts: two low-ma ss X-ray binaries (LMXBs) - AC211 
u;t,e ^L^.(ji9S4); T ,^j l u ia J^Jj4 :1 5 X-2 (jWhite fc Angeliml 
2uuiU DicbaTi et ai. 2005), a dwarf nova c ataclysmic vari- 
able (CV) - M15 CV1 (|Shara et al.ll2004| ). two suspected 
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dwarf novae CVs whose o ptical counterparts are unclear 
(lHannikainen et af 120051) . and a quiescent LMXB - M15 
X-3 (jHeinke et al.l 12009) . In this work we probe into 
these close binary populations through the use of UV 
photometry and color-magnitude diagrams (CMDs). 

Studies of M15's CMD have previously been used to 
examine many features of the cluster's stellar popula- 
tio n. The blue s traggl er (BS) population was examined 
by lYannv et al.l (|1994l) . the hori zontal branch (HB ) and 
giant branches were analyzed bvlCho fc Led ( 2007). and 
UV studies bv lde Marchi fc Parescd (|1994l 119961 ) probed 
into the extremely blue objects, including w hite dwarfs 
(WDs). Most recently, IDieball et all (|2007f ) performed 
a UV analysis of M15, focusing on the variable stars in 
M15's core as well as performing a basic breakdown of 
the UV-bright stellar populations. This study revealed a 
substantial array of "gap" stars that appear in the region 
of the CMD surrounded by the HB, main sequence (MS), 
and WD cooling sequence. This gap zone had previo usly 
been dubbed the "CV zone" bv lKnigge etH] pool due 
to the likelihood of CVs appearing in this part of the 
CMD. 

Optical CMDs of the central region of M15, though 
relatively deep, emphasize populations above the MS 
turn-off (MSTO). With perhaps a few exceptions, it is 
unlikely that C Vs would be visible so far up on an op- 
tical CMD (e.g., lYannv et alj|1994t Ivan der Marel et al.l 
120021) . Observations in the UV, however, result in a dis- 
tinct CMD morphology that highlights populations often 
not included on optical CMDs. In such diagrams, the gap 
zone is emphasized which can result in a larger list of po- 
tential CVs than that which optical filter combinations 
would yield. UV CMDs are also useful in identifying the 
extremely blue features of the HB and WD sequences 
in cluster s. M15 is a low-metallicity cluster ([Fe/H] 
= —2.26; Harris 1996), which results in a correspond- 
ingly blue HB. In addit ion, M15 is known t o possess an 
extended blue HB tail (Moehle r et al.l [19951) and, as we 
discuss in this paper, possibly harbors a small population 
of extreme-HB (EHB) stars which are hotter than the 
hottest HB stars predicted by canonical stellar evolution 
theories. Furthermore, UV CMDs have the unique abil- 
ity to photometrically distinguish a set of He-rich EHB 
stars known as "blue hook (BHk) stars" which are char- 
acterized as being ap parent EHB stars that are unusually 
faint in the UV (e g.. IWhitnev et al.l ll99l ID'Cruz et all 
2000); IBrown et al.ll2001[ ). 

Another intriguing close-binary population expected 
to be found in the central regions of GCs are helium- 
core white dwarfs (He WDs). He WDs are believed to 
arise in cases where mass loss on the red giant branch 
(RGB) is severe enough that the remaining mass in the 
H shell is too small to ignite core He-burning, result- 
ing in an exposed He-core. They are predicted to lie in 
sequences slightly brighter and redder than the canoni- 
cal CO WD cooling sequences d ue to their lower mass 
and therefore larger radii (e.g., iBenvenuto fc Althausl 
1998). The only cluster in which there is a confirmed 
sequence of He WDs is NGC 6397 (IStrickler et al.ll200l 
and references therein) though it has also been suggested 
that a significant fract i on of the WDs in o> Cen are He 
WDs (jCalamida et al.ll2008D . A small number of He 
WDs are also known in each of several other nearby 
clusters as counterparts to ultracompact X-ray binaries 



(jAnderson et all 119931: IDieball et al.ll2005D and m i llisec- 
ond pulsars (lEdmonds et all l200l UFerraro et ail 120031: 
Sigurdsso n et al.l[2003D: in addition one was found in M4 
bv 10 'Toole et al.l (|2006D as a co mpanion to a subdw arf 
B star and one was identified by IKnigge et al.l ([2008D in 
47 Tuc using UV spectroscopy. 

We performed PSF-fitting photometry on optical, 
near-, and far-UV images of M15's central region ob- 
tained by the Hubble Space Telescope (HST) using the 
Advanced Camera for Surveys (ACS). This data set is a 
unique one: confined neither exclusively in the UV nor 
optical regimes; we analyzed M15 using colors that allow 
us to probe deeply into the optically faint yet UV-bright 
popul ations initially uncovered in Ml 5 bv IDieball et al.l 
(2007). By including an optical filter, we further inves- 
tigate the photometric properties of these populations 
with an additional color and the more well-studied opti- 
cal CMD morphology. In Sj3]we describe the observations 
and photometry. In SJ3] we provide a breakdown of the 
various UV-bright populations. Section [4] contains our 
analysis, with ^4.11 containing the description our WD 
cooling models and analysis of the WD population, i j4.2l 
describing our analysis of the radial distribution of pop- 
ulations believed to be of binary origin, and M4.3ti4. 71 fo- 
cussing on the potential EHB, BHk, He WD, and CV 
populations, respectively. Finally, we discuss the impli- 
cations of our results in fJ5] and summarize our results in 



2. OBSERVATIONS & PHOTOMETRY 

Our data consist of HST archival images taken with 
ACS over several epochs in three filters. The filters 
used were F435W, F220W, and F140LP; the first of 
these is similar to Johnson B, while the F220W and 
F140LP bands are in the near- and far-UV, respec- 
tively. The F435W data consisted of 13 separate frames, 
from GO-10401, taken on the High Resolution Channel 
(HRC) all having exposure lengths of 125 s. The F220W 
and F140LP data, from GO-9792, were taken with the 
HRC and Solar Blind Channel (SBC), respectively. The 
F220W data were taken in eight exposures of 290 s for 
a total exposure time of 2320 s. The F140LP data were 
taken in a manner designed to optimize detection of vari- 
ability and consist of 90 individual exposures that result 
in a total exposure time of 24,800 s. The s e two data sets 
are described in detail by IDieball et al.l (|2007l ) . Hence- 
forth, the three filters will be referred to as B435, NUV220, 
and FUVi4 - 

Despite the high resolution of the raw data, the indi- 
vidual frames were insufficient for resolving most stars in 
and around the cluster core. Thus, to increase the reso- 
lution, the image data were combined using a computer 
program developed by one of us (J. A.). The program 
is similar to drizzle with the pixfrac parameter set to 0, 
except that the transformations from each exposure to 
the master frame are based purely on the locations of the 
bright stars — no image header information contributed. 
The resolution of the resulting images was approximately 
doubled and was sampled with a pixel size of 0.0125" x 
0.0125" for all three filters. The central 10" x 10" of the 
stacked images are shown in Figure [TJ the complete field 
of view is 29" x 26" for the HRC and 39" x 31" for the 
SBC. 

We performed photometry using the software pack- 
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age DAOPHOT II (jStetson et al.|[l990h . The detection 
threshold was determined by trial and error for each im- 
age individually in order to ensure that the majority of 
the sources were detected while still limiting the number 
of spurious detections. The IRAF task xyxymatch was 
used to match individual stars across the filters. A tol- 
erance of 2 pixels was used to match between the B435 
and NUV220 frames, while a 3-pixel tolerance was used 
between NUV220 and FUV140 frames. The difference in 
tolerances proved necessary on account of differing geo- 
metric distortions in fields from the HRC and SBC. Sus- 
pected edge detections, spurious detections, and diffrac- 
tion artifacts that passed through this routine were in- 
spected by eye and rejected when necessary. 

Between 130 and 150 individual stars were selected 
from each image as models for a point-spread function 
(PSF). Stars were chosen based on their visual appear- 
ance and the relative crowding of their neighborhoods. 
The PSF was determined in an iterative fashion with 
a careful effort to distinguish between PSF artifacts and 
close neighbors. The hnal quadratically- varying PSF was 
used with the ALLSTAR PSF-fitting photometry routine 
in order to determine the instrumental magnitudes. 

Magnitudes were calibrated and transformed into the 
STMAG system in a similar manner to that outlined 
in detail in t he H ST ACS Data Handbook and by 
iSirianni et allpOOl . Since ALLSTAR performs PSF fit- 
ting photometry, we first had to calculate an initial offset 
between the magnitudes from a small aperture, chosen to 
be 0.05", and the magnitudes produced by ALLSTAR. 
This calibration was performed with the same set of stars 
used to determine the PSF due to the relative brightness 
and isolated nature of these stars. The offset was quite 
uniform across the field, with a standard deviation for the 
calibration stars not exceeding 0.016 magnitudes. Since 
the field is extremely crowded (see Fig. [1} it was impos- 
sible to obtain a magnitude offset for the suggested 0.5" 
aperture. Thus, in order to make the final calibrations 
for our magnitud es, we used the encir cled energy frac- 
tions tabulated in ISirianni et al.1 (|2005[ ) and applied the 
STMAG zero point following the process outlined in the 
HST ACS Data Handbook. 

The calibrate d data set wa s direc tly compared to the 
data set used in lDieball et~aH (|2007D (hereafter D07). Of 
the 1913 stars with both NUV220 and FUV140 detections 
included in the data set used by D07 (A. Dieball, private 
communication 2007) we recovered a unique match for 
1813. In addition, we detected NUV220 sources for 37 
FUV140 detections that had no matched NUV220 source 
in the D07 data set. The 100 stars from D07 that we 
did not recover are not surprising; many of these stars 
are very dim and some lay in regions of the image that 
were not included in our master images (such as the oc- 
culting finger in the F140LP frame) because we relied 
solely on the combined frames for source detection. In 
addition, some of the D07 FUV140 sources appeared as 
two separate sources in our catalog, indicating that we 
possibly resolved some stars that appeared blended in 
their master images. Our final data set contained a total 
of 10,728 individual stars detected in at least two filters, 
necessarily including F220W. This consisted of 3052 indi- 
vidual stars detected in both FUVi 40 and NUV220, 2943 
of which were also detected in B435. Of the 1813 stars 
matched to sources from the D07 data set, 1777 were 



detected in B435 as well. The discrepancy between our 
total number of FUV140 sources, which totaled to 3961, 
and that of D07, which totaled to only 2137, seems to be 
due to our use of a lower threshold for source detection 
which resulted in the inclusion of on the order of 2000 
faint sources that were not included in the D07 data set. 
We expect that most of the false detections in this en- 
larged number of FUV140 sources were filtered out by 
object matching across filters and our visual inspections 
of individual objects. 

While most of our FUV140 and NUV220 photometry 
is consistent with the similar results presented in D07, 
we note that for many of the variable stars our FUV140 
magnitudes are systematically dim as compared to those 
of D07. This is likely a consequence of having performed 
photometry on combined frames that had gone through 
a rejection routine designed to remove bad pixels, warm 
pixels, and other artifacts. This routine rejected any 
pixel inconsistent at the 3 sigma level and therefore would 
affect any stars with variability that exceeded this thresh- 
old. 

3. THE CMDS 

Photometric results are illustrated in the CMDs pre- 
sented in Figures 2-4. In each figure we have identi- 
fied various stellar populations based on their position 
in a given CMD (i.e. FUV140-NUV220 vs. FUV140) and 
then plotted those same stars on a different CMD (i.e. 
NUV220-B435 vs. B435) using a scheme that allows the 
stars, or groups of stars, to be tracked between the dia- 
grams. This technique allows us to better evaluate the 
photometric nature of a star or population by examining 
its position on multiple CMDs and additionally is useful 
for understanding the unfamiliar morphology of the UV 
CMDs. 

We have included multiple versions of the same CMD, 
with different identification schemes for the stellar pop- 
ulations, in order to clearly illustrate the photometric 
properties of each population we analyzed. Individual 
components of the various diagrams are discussed in the 
following sections. Unless otherwise stated, the differ- 
entiation between stellar populations was based on the 
overall appearance of the CMD and the location of the 
stars on the CMD. We caution that this non-quantitative 
method will result in a small amount of confusion be- 
tween the individual groups, but should not have a sig- 
nificant effect on our results. All magnitudes are given 
in the STMAG system. 

3.1. Photometric results for Ml 5 X- 2, CV1, AC211 and 

H05X-B 

We have obtained photometric measurements for 4 
X-ray sources previously identified with close binary 
systems (labeled and plotted as squares in Fig. 2- 
4). These systems are include the confirmed dwarf 
nova (DN) M15 CV1 (hereafter CV1), identified by 
iShara et al.l ( [2004 1 as well as optical counterparts to 
the two LMXBs: AC211, the optical counterp art to 4U 
2127+119 as identified bv lAuriere et al.l (|1984| ). and M15 
X-2 (hereafter X-2), identified as the optical counterpart 
of an ultracompact X-ray binary co rresponding to X- 
ray so urce CXO J212958 1+121 002 bv lWhite fc Angelinil 
(|2001ft and lDieball et al.l (|2005[ ). We also have provided 
photometry in NUV220 and B435 for a faint X-ray source 
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identified in iHannikain en et aT] (|2005l ) as a possible DN 
or a quiescent soft-X-ray transient (qSXT; their "source 
B"). We refer to this source, hereafter, as H05X-B and 
discuss it in detail in §4.5.41 H05X-B was too dim to be 
detected in FUV140 (see Figured]), thus it is not included 
on CMDs which include this filter. 

Photometric data for these four objects are listed in 
Table [2 Our NUV220 and FUV M0 magnitudes for X-2, 
AC211, and H05X-B are consistent with those of D07, 
within the errors. We have also provided photometry in 
B, thus adding to the photometric information for these 
sources. 

CV1 was not detected in the D07 FUV140 master image 
as they found this star too faint and too near a bright 
neighbor to be detected. But, they do detect the star 
during an outburst phase in one of the observing epochs, 
therefore confirming it as a DN. In our FUV140 mas- 
ter image, the source corresponding to CV1 is difficult 
to resolve by eye but the FIND task in DAOPHOT II 
was able to identify it as a discrete source (see Figure 
[6]). However, when run through the PSF-fitting photom- 
etry of ALLSTAR, the proximity to the bright neighbor 
proved an issue and caused the star to be subtracted as 
part of the PSF of the neighboring star; thus we used 
a FUV140 magnitude derived from aperture photometry. 
We report the NUV220 and B435 photometry of CV1 with 
equal confidence to stars of similar magnitudes, but note 
that since we were not able to use PSF-fitting for our 
FUV140 photometry of CV1 the detection may be more 
questionable and our results less reliable. 

The NUV220 magnitude for CV1 reported by D07 is 
significantly brighter than our reported magnitude. As 
noted in the erratum to that paper, this is likely related 
to the fact that we performed PSF- fitt ing photome- 
try w hile D07 used aperture photometry (jDieball et al.l 
120101 ). With aperture photometry, a significant amount 
of flux from the bright neighbor star would likely end 
up in the aperture for the CV1 counterpart, and thus 
cause the calculated magnitude to be too bright. When 
we performed aperture photometry with an aperture of 
comparable size, it resulted in a magnitude consistent 
with that in D07. 

3.2. The Main Sequence and Red Giant Branch 

The MS extends approximately 3.5 magnitudes below 
the turn-off in B435 (see small orange dots in Figure [2£i). 
The MS has a width of about 0.5 magnitudes just be- 
low the turn-off and widens somewhat at the faint end. 
There is scatter to either side of the MS, some of which 
likely is due to photometric error, but may also indicate 
potential binary systems. The detection limit at B435 ~ 
22.5 and NUV220 - B435 ~ 2 is due to the depth of the 
NUV220 data. Since MS stars are relatively red, there 
are many fewer FUV140 detections of MS stars and the 
MS in Figure Oi extends only about one magnitude be- 
low the turn-off with a significant amount of scatter due 
its proximity to the detection limit. 

The RGB (large red dots) is very well defined in Fig- 
ure [2ji, forming a tight, narrow sequence stretching over 
4.5 magnitudes in B435. The FUV140 filter strongly sup- 
presses the presence of the red giant population and the 
fact that a large number of RGB stars are detectable in 
the FUV140 frame at all is a bit of a mystery. It is most 
likely due to the "red leak" phenomenon associated with 



the U V fil ters and discussed in iChiagberge fc Siriannil 
pOOl and iBoffi et ail pOOl . Red leak can greatly af- 
fect photometric measurements of red stars, but since the 
focus of this paper is UV-bright stars it should not af- 
fect our results. Because it is very difficult to distinguish 
MSTO stars from RGB stars in the FUVi 40 - NUV220 
vs. FUV140 CMD we make the distinction only using the 
NUV220- B435 vs. B 435 diagram. 

3.3. The Asymptotic Giant Branch 

The asymptotic giant branch (AGB; small teal dots) is 
most easily identified in Figure [5^ as the slightly curved 
sequence turning away from the tip of the RGB and span- 
ning about 2.5 magnitudes in color. The transition be- 
tween the AGB and the HB is somewhat unclear; for the 
purposes of this paper we have defined the edge of the 
AGB just redward of the clump of variable stars belong- 
ing to the RR Lyrae instability strip. Due to this loose 
distinction between the HB and AGB, there may be some 
confusion between the blue edge of the AGB and the red- 
dest HB stars, however this is not of particular concern 
as these stars are not the focus of this paper. 

Like the RGB, the AGB is not as well defined in Figure 
[Sh. because very little flux from these stars is emitted 
in the far-UV bandpasses. However, the AGB can be 
identified as the sequence of stars connecting the tip of 
the RGB to the HB. In Figures Wp and [3]a there is a 
clear bend toward the blue at the point where the AGB 
transitions to the HB (FUVi 40 ~ 22.5). 

3.4. The Horizontal Branch 
3.4.1. Normal Horizontal Branch Stars 

The main portion of the HB (plotted as large light 
blue dots on Fig. 2-4) consists of a small set of red 
horizontal (RHB) stars and much larger set of blue 
horizontal branch (BHB) stars separated by variable 
stars in the RR Lyrae instability strip. A significant 
RHB population is not necessarily expected since metal- 
poor stars produce BHB stars; however, the presence of 
the RHB population i n M15 is well docu mented (e.g., 
iBuonanno etall H985T : IPreston et al.l l2006h . The BHB 
begins at the blue edge of the RR Lyrae strip extend- 
ing toward the WD sequence. For comparison we have 
plotted a theoretical zero-age horizontal branch (ZAHB; 
solid dark red and dark blue lines) generously provided 
by Santi Cassisi ( private communic a tion, 2009 & 2010) 
and described in iPietrinferni et al.l (|2004l ). on Figures 
[2] - 0J They were generated assuming a metallicity of 
Z=0.0001 (corresponding to an [Fe/H] —2.3) and were 
ad justed for exti n ction following the method outlined 
in iCardelli et al.1 (11989T). We assumed a color excess 
of E(B-V)=0.1 (|Harris I fl996h which resulted in a cor- 
rection A\ = 0.81±0.11, 0.94±0.12, and 0.42±0.08 for 
FUV140, NUV220, and B435, respectively . 

For the dark blue curve , a distance of 10.3±0.4 kpc 
(|van den Bosch et al.l 120061) was assumed resulting in a 
distance modulus of 15.06. This is in general the most 
accepted distance measurement for Ml 5, however, it is 
clear in the CMDs including the FUV140— NUV220 color 
that using this distance modulus does not produce a good 
fit for the red side of the HB, as the ZAHB appears sig- 
nificantly brighter than the observed HB stars. If the dis- 
tance modulus is adjusted to 15.25, consistent with the 
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distance modulus determined by Kraft fc Ivans! (|2003[) . 
(solid dark red curve), a much better fit is achieved. 
However, with the larger distance modulus the ZAHB 
is significantly dimmer than expected on the blue side of 
the HB. So, neither distance modulus produces an ideal 
fit for the entire HB and we ascribe this to the limita- 
tions of these models in fitting our dataset. Yet, since 
the fit is reasonable in both cases, we believe the mod- 
els to be sufficient for the purposes of this paper. The 
dotted lines associated with each curve represent the la 
error bars including both the erro r in distance and red - 
dening. For the dark blue curve, IKraft fc Ivaiisl (2003) 
do not cite an error, so we assumed a reasonable error of 
0.10 magnitudes in the distance modulus. 

3.4.2. Extreme Horizontal Branch Stars 

As previously noted b y many authors 
Ce.g.. iDorman et all 119931 iFerraro et all 119981 : 
iDalessandro et all 120081) the hottest BHB stars owe 
their extremely blue colors to thin H envelopes which 
result in higher effective temperatures. The thin enve- 
lope is believed to be due to significant mass loss on 
the RGB, likely to a close binary companion or possibly 
stellar winds. Because of the thin envelope these stars 
do not undergo the usual evolutionary progression to 
the tip of the AGB. After leaving the HB, they have an 
insufficient mass in their hydrogen shell to enable the 
onset of the thermal pulsation phase; instead they either 
move off the AGB early and become post-early AGB 
(P-EAGB) or never ascend the AGB at all, becoming so 
called AGB-manque stars. 

This group, composed of the hottest BHB stars, is 
known as the extreme horizontal branch (EHB) and is 
often defined as the set of BHB stars with effective tem- 
perature greater than 20,000K (e.g.. IBrown et al.l 120011 : 
IMoehler et aHl200i IDieball et al] l2009) but is more dif- 
ficult to define observationally. Spectroscopically EHB 
stars correspond t o subdwarf B (sdB) and subdwarf OB 
(sdOB) field stars (Hcbcr 19871). but photometric defini- 
tions differ between authors. For clarity, we use the pre- 
viously mentioned ZAHB models, and chose stars with 
T e s > 20,000K as EHB stars; using this method we iden- 
tify 6 stars as EHB candidates (maroon pinched squares 
on Fig. 2-4). One of these candidates was not detected 
in B435 due to its location near the edge of the B435 im- 
age, so all 6 candidates only appear on Figures [2p, [3^, 
and|4p and only Figure [3^ displays all 6 candidates with 
the expected shape and color. 

3.4.3. Blue Hook Stars 

Five of the 6 potential EHB stars appear subluminous 
in the UV compared to BHB stars with similar temper- 
atures (plotted as EHB stars in diamonds in Fig. 2-4). 
This indicates these as candidate blue hook (BHk) stars. 
(The subluminous nature of these 5 stars is perhaps most 
easily seen in Fig. 2}d.) Defining EHB and BHk stars 
photometrically can be a very difficult task as there is 
no clear consistent photometric definition. Here we have 
chosen to rigorously define EHB stars as only those with 
Tcff > 20,000K and have chosen our BHk candidates as 
those that appear dimmer than the ZAHB models of sim- 
ilar temperature HB stars in the CMDs containing the 
FUV140— NUV220 color. Clearly this definition depends 
somewhat on which ZAHB curves we use. As can be 



seen in Figures [5p, |3K, and Hp, four of the BHk candi- 
dates lie fairly clearly below the canonical ZAHB while 
the last one appears subluminous compared to the ZAHB 
using the smaller distance modulus (dark blue curve) but 
is more well within the error bars of the dimmer curve 
which uses larger distance modulus (dark red curve). We 
still consider it a possible BHk star but recognize that its 
identification as such is weaker than the other four. 

Data on the stars we have identified as BHk candidates 
can be found in Table [21 which also contains photometric 
data for other UV bright sources whose classification is 
uncertain. These stars, as well as all the Table 2 stars, 
are plotted on Figure [7J 

BHk stars are similar to EHB stars but have helium 
rich envelopes that cause increased opacity in the at- 
mosphere and leads to their subluminous nature. They 
are an intrinsically rare and relatively poorly under- 
stood populatio n that have only recently begun to be 
identified (e.g.. IBrown et al.l 120011: IBusso et all 120071: 
iSandauist fc Hessl 120081 ). At least two scenarios have 
been presented as to the origi n of BHk s tars. One sce- 
nario, originally proposed bv iLee et all (|2005l l is that 
BHk stars and the hottest HB stars may result from 
the standard evo lution of a He-rich subpopulation in the 
cluster (see also |Dieball et al.l 120091 : iD'Antona fc Calol 
2008;; IMoehler et alll2007D. The other scenario discussed 
by many authors (e g.. IBrown et al.l 120011: Ib'Cruz et all 
120001: IMoehler et al.ll2004D suggests that BHk stars lose 
a significant amount of mass along the RGB, likely to 
a close binary companion, such that they do not ignite 
helium at the tip of the RGB, but instead, undergo a 
helium-core flash as they begin to descend the white 
dwarf cooling sequence - a so called late He-flash. 

The "blue hook" terminology is a reference to the ap- 
pearance of this population on UV CMDs. While M15 
does not currently have a recognized population of BHk 
stars, one helium -rich hot BHB star was spectroscopi- 
cally identified bv IMoehler et al.l (|1997l) that would most 
likely appear as a BHk star on a CMD. This star is well 
beyond our field of view, so we can provide no further 
insight on it. However, the 5 candidates we have iden- 
tified could establish the presence of a BHk population 
in Ml 5, and may provide an important clue about the 
origin of these stars. The significance of BHk stars in 
M15 is further discussed in gT3]& El 

3.5. White Dwarfs 

A prominent population of WD candidates that is 
in general agreement with the theoretical cooling mod- 
els is apparent in the CMDs in Figures 2-4 (magenta 
pinched triangles). We have detected 64 WD candidates 
in FUV140 and only 25 in B435 (20 of which were also de- 
tected in FUV140). The majority of the WD candidates 
detected in FUV140 were not detected in B435 due to the 
faintness of these stars in B435 and crowding problems 
that become particularly pronounced in the B435 image. 
As can be seen by examining the cooling curves plot- 
ted on Figures 2-4, some WD candidates appear to be 
more consistent with the cooling sequences for helium- 
core WDs (He WDs) than for normal carbon-oxygen 
WDs (CO WDs). This, along with the details of the 
WD population and models, is discussed in §4.11 

3.5.1. Helium-Core White Dwarfs 
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He WDs are the remnant of a star that loses most of 
its hydrogen envelope on the RGB and therefore never 
undergoes a He-flash or subsequent HB and AGB evolu- 
tion, but instead ends up cooling as a degenerate helium 
core with a thin hydrogen envelope. He WDs are usually 
found as members of a binary system and are believed 
to arise primarily due to Roche lobe overflow and m ass- 
transfer in a close binary system (e.g., Wcbbink 1975). In 
dense environments such as GC cores it is possible that 
He WDs are formed by collisions involving RGB stars 
that result in a commo n envelope phase a nd eventual 
ejection of the envelope (IDayies et al.l 119911 ) . Although 
iCastellani fc Castellanil (|1993l ) show that in cases of ex- 
treme mass loss from winds it is possible for isolated stars 
to become He WDs, this seems an unlikely scenario for 
the production of a substantial number of He WD as the 
amount of mass loss required is much larger than can be 
explained by canonical stellar evolut ion alone. 

Based on theoretical studies, (e.g.. [Hansen et aLll2003f) 
GCs, especially those that have undergone core col- 
lapse, should possess observable He WD sequences. How- 
ever, only one GC (NGC 6397) has an extended He 
WD sequence c urrentl y identified (e.g..lCool et al.lll998t 
Edmonds et al.l I1999{: iTaylor et al.l 120011: Hansen et al.l 
2003t iStrickler et all 



been previously identified in M15, but in §4.1l we present 
strong evidence for the presence of a He WD sequence 
and further discuss this population in §4.41 

3.6. Blue Stragglers 

The BS sequence appears as an extension of the MS to 
luminosities greater than the turn-off point. It is gen- 
erally believed that these stars are the products of a 
merger of two or more MS stars that, upon merging, pro- 
duce a core hydrogen-burning star more massive than 
the MSTO. There is an ongoing discussion of whether 
the most important formation mechanism for BSs is 
the gradual coalescence of a close bina ry systems or 



physi cal collisions between two stars (e.g.. [Fregeau et al 
Bo llMapelli et al.ll2006t iLeigh et al.ll2007tlKnigge et al 



20091 ). However, it is well agreed that the BSs found pref- 
erentially in the dense central regions of clusters. 

We have identified 31 BS candidates from their position 
in Figure [2ji (small blue inverted triangles). However, 
it is difficult to judge where the MS ends and the BS 
sequence begins. It appears that Figure [3^ may be more 
insightful for distinguishing BSs from MSTO stars as the 
sequence of stars bluer than the MSTO is stretched out, 
thus making BSs more clearly distinct from the turn-off 
stars. Using Figure |3^, we have identified approximately 
53 BS candidates. 

3.7. Gap Objects 

We also have identified a significant set of stars that we 
term gap objects following D07 (plotted as bright green 
dots and open circles in Figures 2-4) . This population is 
composed of stars that populate the gap zone between 
the MS and WD regions of the CMD where there is an 
increased likelihood of finding CVs. The previously iden- 
tified cataclysmic variable, CV1, appears in this region 
as expected. 

In Figure [3£i we identified 60 stars clearly populating 
the gap zone (thus classified as gap objects) however in 



Figure^ we identify only 22 gap objects, most of which 
lie near the MS or the faint-end photometric limits. Most 
of the sources identified as gap objects in Figure [3£i were 
detected in B435 but rather than appearing in the gap 
region on NUV220 - B435 CMDs they appear primarily 
on the MS. This puzzling feature is discussed in detail 
in £14.51 . We are confident in the significance of the gap 
population as identified in Figure [3^, as these sources 
were all inspected by eye and D07 found a very similar 
population in their FUV140 and NUV220 photometry. 

3.8. Bright Blue Gap Objects 

The population of stars that we identified as bright 
blue gap objects (pink asterisks in Fig. 2-4) consists of 
8 stars (including X-2) detected in each filter that are 
bluer than the MS in both NUV220- B 435 and FUVi 40 - 
NUV220 by at least 1 magnitude but are found be- 
tween the standard WD and HB sequences in bright- 
ness. Canonical stellar evolution does not include stages 
expected to populate this region of the CMD for any 
significant time. WDs rapidly transiting from the post- 
AGB phase toward the WD cooling sequence may pass 
through this region, but the timescale is such that the 
likelihood of detecting even one star in such a phase is 
extremely small. Thus, from a canonical stellar evolution 
standpoint, the presence of several objects in this region 
is unexpected. In Figures 2-4 it is however clear that 
the theoretical He WD cooling sequences run through 
this region and we will discuss the plausibility that these 
stars may be young He WDs in §4.4.41 

This population is perhaps most clearly defined in Fig- 
ure GJi as the stars with with FUV140— NUV220 < —0.5 
and 16.4 < FUV140 < 18.4. Using this criteria we have 
identified 10 bright blue gap objects; however, in Fig- 
ure [3Jd, one appears more convincingly to be a hot WD 
and another appears on the RGB at NUV220 - B435 
4. The nature of the latter star is quite confusing; a 
plausible explanation is that we have detected a chance 
superposition of an RGB star and a very blue object, 
such as an HB star. When examined by eye, the cen- 
ter of light appears to be slightly inconsistent across the 
three images, supporting this hypothesis. 

It is unclear what the true nature of these bright blue 
gap objects is and we consider several possibilities in §33 
One of the 8 objects that appear as bright blue gap ob- 
jects in both colors is X-2, raising the possibility that 
some of these stars could be close binaries with accretion 
disks. These stars may also be related to the BHk or 
He WD populations and each of these possibilities are 
discussed in i j4.7l Stars identified as bright blue gap ob- 
jects are included in Table [5] and Figure [7] as UV8 and 
UV10-12. 

4. ANALYSIS 

4.1. White Dwarf Cooling Curves 

A set of cooling sequence models for both CO WDs 
(solid blue curves) and He WDs (dashed green curves 
and dotted purple curves) have been calculated and plot- 
ted on Figures [8] & M Figures 2-4 also include the CO 
WD cooling sequences (solid lines) and one set of the 
He WD cooling sequences (dashed lines) for orientation 
purposes. The code and inp ut physics for t hese mod- 
els are described in detail by iSerenelli et al.l (|2002D and 
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are calculated in the same manner as the models used in 
the analysis of the H e WD population in NGC 6397 by 
iStrickler et all (|2009f) . The models cover a mass range 
of M = 0.45 - 1.10 M for the CO WDs and M = 0.175 
- 0.45 M Q for the He WDs (masses indicated in caption 
for Fig. EI&EJ). The models in Figures M & M differ from 
each other in that the He WD models in Figure 8 (purple 
dotted lines) have thin hydrogen envelopes and the He 
WD models in Figure 9 (green dashed lines) have thick 
hydrogen envelopes. The thickness of the hydrogen en- 
velope affec ts the color and cooling times for He WDs a s 
discussed in lSerenelli et all (|2002f ). lAlthaus et al.l (|2001f ). 
and briefly in §4.4 | of this p aper. A distance of 10.3 kpc 
(|van den Bosch et al.ll2006l) was assumed and the cooling 
curves were adjusted for extinction in the same manner 
as the ZAHB. We have included error bars in the upper 
portion of Figures [8] & [9] that represent the la error in 
both distance and reddening for the cooling curves. 

The He WD models have a progenitor metallicity of 
Z=0.0002 (corresponding to [Fe/H] ~ —2), consistent 
with the derived metallicity for M15. The CO WD mod- 
els are all generated for progenitor stars with solar metal- 
licity. Since the evolution of CO WDs depends very lit- 
tle on nuclear burning and the cooling timescale is not 
metallicity dependent, we consider these tracks reason- 
ably suitable for M15 despite the metallicity discrepancy. 

From their location on the CMDs in Figures 2-4, we 
have identified a total of 73 stars that appear WD-like 
in at least one CMD. These stars are plotted as dots and 
open circles on Figures [8] & GO It is somewhat difficult 
to distinguish CO and He WD candidates as there is 
large uncertainty in the photometry for the dim stars. 
Table [3] contains a summary of our assessment of these 
WD candidates as discussed in detail in the following 
paragraphs. 

Using their location on the CMDs, we determined 11 
of these stars represent likely candidates for CO WDs, 45 
appear to be likely He WD candidates, 10 appear to be 
good WD candidates but are ambiguous as to whether 
they belong to the CO or He WD population, and 2 
are variable stars from D07 (plotted as filled triangles in 
Fig. [8] & [9]) that appear on the CO WD curves, but are 
more likely CV candidates. The remaining 5 represent 
possible WD candidates but are subject to significant 
photometric scatter that makes their nature unclear. 

Of these 73 stars, only 20 were detected in all three 
frames (see filled black and red dots in Fig. |8]&[9]). For 
3 of these it is unclear whether they are truly WDs: they 
appear WD like in Figures [5b & [3b but lie more than 
0.5 magnitudes to the red side of the tracks in Figures 
[8ji & [9^, thus they may actually be CVs. This leaves us 
with 17 stars which we consider to be our strongest WD 
candidates because we have two-color data that allows 
us to better distinguish their true nature. 

We draw the following conclusions concerning these 17 
stars: (1) Five are likely candidates for CO WDs. Two 
are strong candidates that appear on the CO WD cooling 
curves in both diagrams and the other three are consis- 
tent with being either low-mass CO WDs or young, mas- 
sive He WDs. (2) Seven are good candidates for He WDs 
because they appear reasonably consistent with the He 
WD cooling sequences in both diagrams and are clearly 
separated from the CO WD cooling curves. These are 
the 7 stars plotted as large red filled circles on Figures 



IH&IHl (3) Five of the 17 stars are clearly WD-like but 
are ambiguous as to which population they belong since 
they appear on different cooling curves in each diagram. 

We caution that, for many of these stars, the distinc- 
tion between He WDs and CO WDs is very difficult and 
subject to assumptions such as the cluster reddening and 
distance (especially in CMDs using FUV140— NUV220 as 
these two quantities are the most sensitive to redden- 
ing). The effect of the uncertainties in these quantities 
can be seen by examining the la error bars for the cooling 
curves which are plotted in the upper portion of Figures 
[8] & [9l If we assume less reddening, some He WD can- 
didates become better CO WD candidates, and if we as- 
sume more reddening, some CO WD candidates become 
He WD candidates. Despite these uncertainties, we still 
find our claim that M15 possesses some population of He 
WDs to be strong. 

In addition to the WD candidates discussed above, the 
bright blue gap objects are plotted on Figures|S]&in](grey 
asterisks) because a subset of these objects appear to 
be consistent with the thick H envelope He WD cooling 
curves. Of the 7 bright blue gap objects detected in all 
three frames (not including X-2), all appear to be con- 
sistent with being young, low-mass He WDs with thick 
hydrogen envelopes (see Fig. |HJ). The nature and origin 
of He WDs, the significance of the thin versus thick H 
envelope models, and the implications of the existence of 
He WDs in M15 is discussed in more detail in £14.41 

4.2. Radial Distributions 

Binary systems with a total mass larger than the av- 
erage stellar mass in the cluster are expected to segre- 
gate towards the cluster core on the time scale of a half- 
mass relaxation time due to dynamical friction. Since 
M15 is a core-collapsed clu ster, with a h alf-mass relax- 
ation time of i r ,hm ~ 1 Gyr (Harris 1996), mass segrega- 
tion should have taken place and any "massive" binaries 
should be centrally concentrated. Many of the popula- 
tions addressed in this paper (BS, BHk, EHB, He WD, 
CV) are likely members of such binary systems, thus are 
expected to be segregated towards the cluster center. 

In Figure [10] we have plotted the cumulative radial dis- 
tributions for the BSs, gap objects (which have been sep- 
arated into two subsets, see ^4.2. 3p . and bright blue gap 
objects, along with the radial distribution of HB stars 
for comparison. We chose the HB population as the ref- 
erence population because it is bright enough that com- 
pleteness issues should be minimized, even in the dense 
central regions. To determine the statistical center of the 
cluster we used the MSTO population from Figure^ us- 
ing the positions in the NUV220 image. We must note 
that any dim population will be drastically incomplete 
in the most central part of the cluster of the FUV140 
image because the extended PSF haloes and diffraction 
artifacts from several bright stars concentrated near the 
cluster center effectively mask many dim stars in this 
region (see Fig. [lj. 

In order to analyze the significance of any apparent 
central concentrations, we performed Kologorov-Smirnov 
(KS) tests comparing the radial distribution of each pop- 
ulation to our reference HB population. The KS statistic 
gives the probability that the two samples being com- 
pared are consistent with being drawn from the same dis- 
tribution; therefore, a lower KS probability corresponds 
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to a more statistically significant difference between the 
distributions of the two samples, thus a more centrally 
concentrated sample. 

4.2.1. Distribution of Bright Blue Gap Objects 

We hnd a KS probability of 0.04% and 0.4% for the 
bright blue gap objects (as identified in Figure [2^ &[3K, 
respectively) as being consistent with the HB sample. 
This is a strong indication that these stars are more mas- 
sive than the average HB mass in the cluster, likely due 
to a binary nature. 

4.2.2. Distribution of Blue Stragglers 

We find a KS probability of 1% that the BS population, 
as identified from Figure [5^, is consistent with the HB 
sample. This indicates a significant segregation towards 
the cluster center. Yet, for the BSs identified from Figure 
[2^, the KS probability is 30%. This is a surprising result 
considering that the latter set represents the brightest 
and presumably most massive BSs, which are expected 
to display the strongest effects of mass segregation. Also, 
nearly every BS in Figure [2^l also was chosen as a BS in 
Figure |3K, implying that it is the dimmer, less massive 
BSs that show the strongest degree of segregation; upon 
further inspection we find this to be true. 

If we consider the stars identified as BSs in Figure 
that were detected in all three frames and divide them 
into a "bright" group with B 435 < 18 and a "dim" group 
with B 435 > 18, we find that the "bright BSs" have a KS 
probability of 13% while the "dim BSs" have a KS proba- 
bility of 1%. This would seem to indicate that these dim 
BSs are significantly more centrally concentrated than 
the bright BSs. Although a rigorous explanation for this 
phenomenon is beyond the scope of this paper, we sug- 
gest it may be a statistical issue due to the smaller num- 
ber of BSs in Figure [2] (only 31 BSs selected from Fig. [2] 
vs. 53 selected from Fig. [3]). 

Recently, a significant number of globular clusters have 
been discovered to have a bimodal BS distribution with 
a peak inside the core radius then a quick drop off in 
a region termed the "zone of avoidance " followed by 
a second peak at several core radii (e.g.. iFerrarq et al.l 
H99l IMapelli et all 12005 iDalessandro et al.l 120081 ). We 
don't expect to see such a distribution in our data since 
the expected radius for the drop off in M15 is well be- 
yond our field of vie w. Following the method outlined in 
IMapelli et al.l (|2006[ ) with reasonable parameters for M15 
we estimate a zone of avoidance at r ss 3.3 — 7.5' (10-20 
pc), depending on the exact parameters used in the cal- 
culation. Even though this is a relatively large range it 
is undoubtably beyond our field of view and near the 
half-mass radius of the cluster. 

4.2.3. Distribution of Gap Objects 

For the gap objects we find a KS probability of 89% 
(from the Figured group) and 80% (from the Figure [3^ 
group) as being consistent with the HB. This indicates 
that the gap objects do not show a central concentration. 
But, since many of the gap objects are faint, this popu- 
lation suffers serious incompleteness issues, especially in 
the innermost regions of the cluster. In an attempt to 
lessen this bias we have plotted the distribution of only 
those gap objects with NUV220 < 22. The gap objects 



are plotted as two subsets on Figures 2-4 (NUV220 < 
22: large bright green dots; NUV220 > 22: bright green 
open circles), and thus it can be seen that by using this 
criteria alone we were able to select a group that should 
have completeness levels that are more consistent with 
the other populations we are analyzing. 

For this brighter subset (NUV220 < 22) we find a KS 
probability of 8% and 64% from the Figure [2] and [3l re- 
spectively. These are still relatively high KS probabilities 
and therefore do not necessarily indicate the presence of 
a statistically significant central concentration for these 
populations, but they do suggest that the brighter set of 
the population derived from Figure ®l may be centrally 
concentrated to some degree. Despite the moderate KS 
probabilities we report, we still consider it likely that at 
least some of these objects are close binary systems based 
on their photometric properties and suggest that the KS 
probabilities given here be regarded as upper limits due 
the completeness and crowding issues discussed above. 

4.2.4. Distribution of Other Populations 

Both the EHB and He WD populations are possibly 
members of binary systems more massive than the aver- 
age stellar mass in the cluster and therefore may be seg- 
regated to the center as well. We, however, do not report 
on the radial distribution statistics for these populations 
due to issues of small number statistics and completeness, 
respectively. The EHB is composed of only 6 stars, so 
the results from any statistical analysis we perform would 
not be robust. The He WD population has enough can- 
didate members to allow for more robust statistics, but 
the dim nature of these stars makes them suffer drasti- 
cally from the completeness and crowding issues in the 
cluster center that have been discussed throughout this 
paper. Therefore, there is an apparent dearth of these 
stars in the innermost region of the cluster which is al- 
most certainly due to crowding issues that mask these 
dim stars. However, because of this feature, our data 
is not sufficient to provide meaningful statistics on their 
radial distribution. 

4.2.5. Comparison with D07 

D07 also performed an analysis of the radial distribu- 
tion of the BSs, and gap objects compared with the HB 
stars. Consistent with our results, they found that the 
BSs are significantly more centrally concentrated than 
the HB. However they find a KS probability of the two 
samples being from the same parent distribution of 14%, 
which is much larger than the KS probability that we 
report of 1%. This is most likely due to differences be- 
tween the stars classified as BSs; D07 classified 69 stars 
as BSs while we only classified 53 stars as such. 

Our finding that the gap objects do not necessarily 
appear to be centrally concentrated seems at odds with 
the finding in D07 that the gap population is the most 
centrally concentrated with a KS probability of only 8% 
when compared with the HB stars. However, D07 only 
considered stars with NUV220 < 21, so their result should 
only be compared with our result for the "bright" gap ob- 
jects with NUV220 < 22. While the samples are still not 
entirely similar, it is striking that we find a much larger 
KS probability of 64%, even for this brighter sample. 
This could be, in part, due to the fact that we extend 
our cut a magnitude deeper and therefore include more 



UV-Bright Stellar Populations in M15 



9 



stars and may have more issues concerning completeness 
near the cluster center. However, a large contribution to 
this difference is almost certainly due to the fact that the 
gap population in D07 included the stars we have termed 
bright blue gap objects which are very strongly centrally 
concentrated. 

4.3. Blue Hook Candidates 

Of the many UV-bright stellar exotica that have been 
discussed thus far, evidence of a BHk population is of 
particular interest as it may provide an important clue 
about how these stars originate. The presence of BHk 
stars in Ml 5 is somewhat contrary to the hypothesis 
that BHk stars are the product of a He-rich subpop- 
ulation, because with a m ass of only 4.4 x 10 5 M 
(|van den Bosch et al.l [2006) it is unclear whether M15 
has a potential well deep enough to maintain sufficient 
gas from stellar ejecta to form a second generation of 
stars. In addition, M15 does not show a split MS as ex- 
pected for a cluster containing a second generation (e.g., 
NGC2808: Piotto et al. 2 007; u Cen: Bedin et al. 2004). 
ID'Antona fc Caloil (|200l argue that the period distribu- 
tion of RR Lyrae stars in M15 is best explained by the 
existence of a He-rich subpopulation, but note that their 
models produce HR diagrams which are not particularly 
good fits for the morphology of M15's HB. So, while a 
subpopulation in M15 should not be ruled out, it seems 
a somewhat unlikely explanation for BHk stars in this 
cluster. 

Further more, the high central densi ty of M15, 7 x 10 6 
M Q pc -3 (jvan den Bosch et al.1 120061 ) . results in a very 
a high interaction rate weighing in favor of the late He- 
flash origin. In this picture the He-flash occurs on the 
WD cooling sequence, thus it is necessary that there be 
enough mass loss during RGB phases to avoid the He- 
flash at the tip of the RGB . The high central density and 
high interaction rate of M15 should lead to an increased 
number of mass-transfer binaries and collisions that may 
strip the outer layers of an RGB star, thus increasing the 
potential for BHk stars to be formed this way. 

4.4. Helium White Dwarfs 

We have identified a substantial number of He WD 
candidates (see §4.1j) . However many of our candidates 
lie very near our detection limit and therefore suffer from 
significant photometric scatter and are more susceptible 
to being false detections. For this reason, we will focus 
our analysis on those 7 He WD candidates that were 
detected in all three filters and appear as good candidates 
for He WDs in both Figures |8] & |9] (plotted as large red 
filled circles). 

4.4.1. Thick vs. Thin H Envelopes 

One question of current interest about He WDs is how 
massive their hydrogen envelope is. The candidates iden- 
tified here, in conjunction with the models discussed 
in Sj4Tj may provide some insight on this parameter. 
The mass, and thus thickness, of the hydrogen envelope 
greatly affects the cooling timescale and, to a lesser ex- 
tent, the photometric properties of a He WD of a given 
mass. Since it is unclear what envelope mass is expected, 
we have included models for both "thick" and "thin" en- 
velopes. In the thick H envelope case nuclear burning 



(pp-chain) at the base of the envelope dominates the ra- 
diation and cooling time, while in the thin H envelope 
case the contribution from nuclear burning is negligible 
so the remnant cools via thermal radiation. The mod- 
els presented here represent the two opposite extremes in 
the the range of hydrogen envelope ma ss and therefore, 
as pointed out by IStrickler et al.l (|2009D . should "bracket 
reality." For further details of the models and the im- 
plications for "thick" vs. " t hin" e nvelo pes the reader 
is refe rred to ISerenelli et al.l (|2002j ) and iStrickler et al.l 
(|20M) . 

4.4.2. Photometric Masses 

It is difficult to derive exact photometric masses from 
our candidate WD populations, however in Figure [9] the 
7 strong He WD candidates seem to be well bracketed by 
the curves spanning the mass range from M « 0.200 M - 
0.275 Mq. The 7 bright blue gap objects (grey asterisks) 
that appear to be He WD candidates in Figure |H] also 
appear to be most consistent with the models in this 
mass range. The photometric mass range is much less 
well constrained in Figure HI as the set of 7 strong He 
WD candidates seem to span the mass range from M w 
0.175 M - 0.350 Mq. 

It is nearly impossible to approximate a photomet- 
ric mass using the entire set of 45 He WD candidates. 
The candidates span the entire range of model He WD 
masses in Figure |H] and approximately 10 of the candi- 
dates appear significantly redder than the reddest model 
in Figure [H It is possible that the reddest stars are not 
truly He WDs, but instead are cither CVs or detached 
WD-MS binary systems that just happen to lie near the 
WD cooling sequence. However, since this region is rela- 
tively well populated and reasonably well separated from 
the gap object population in color, we still present them 
as more likely WD candidates than gap objects. This 
still does not eliminate the possibility that some of our 
candidates may be CVs, detached WD-MS binaries, or 
even LMXBs especially as this region contains two vari- 
able stars classified as likely CVs based on the variability 
study by D07. 

4.4.3. Cooling Ages and Implication for Formation 

Following the analysis in IStrickler et al.l (|2009l ) we an- 
alyze the apparent cooling ages of the He WDs in order 
to gain insight into the possible formation mechanisms as 
well as investigate the plausibility of the thin vs. thick H 
envelope models. Focussing on the 7 strongest He WD 
candidates in the thin H envelope case (Figure [5]), we 
see that all 7 stars have cooling age <100 Myr, with the 
dimmest having a cooling age approximately equal to 100 
Myr. Therefore we can calculate an average formation 
rate of about 70 Gyr -1 . Yet 4 of these stars have cooling 
ages of < 2 Myr which implies an implausible formation 
rate in recent epochs of 2000 Gyr -1 . 

We can estimate the rate at which stars are turning off 
the MS in Ml 5 for comparison to the formation rate im- 
plied by the cooling ages of our WD candidates. Adding 
our sample of MSTO stars that are clearly redward of 
the turn off point (NUV220 - B435 ss 1.5) to our number 
of RGB stars, we have a total of about 800 stars in these 
phases. Given that the time that a star spends in such 
phases is approximately 1.5 Gyr (|Pols et al.l I1998T) . we 
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obtain a rate for stars turning off the MS of 530 Gyr -1 . 
This means that the formation rate implied by the 4 
stars with cooling ages of < 2 Myr using the thin H 
envelope models is significantly larger than the rate at 
which stars are leaving the MS and for this reason alone, 
we find that the formation rate implied from the thin H 
envelopes models to be unreasonable. It is possible that 
these brighter candidates arc not He WDs but are in- 
stead WD-MS detached binaries or CVs, so we can not 
rule out the thin H envelope model altogether with just 
this evidence. 

Additionally, we should note that the estimation of 
the rate at which stars are turning off the MS can be 
calculated in several ways and may be sensitive to issues 
that are difficult to account for, such as mass segrega- 
tion. Another way to calculate this rate is to use the 
HB; using the total number of stars found in the HB (w 
130) and an estimate for the lifetime of a HB star 100 
Mvn lCassisi et al] |2004) we get a turnoff rate of approx- 
imately 1300 stars Gyr -1 . This is about 2.5 times larger 
than the rate obtained using the MSTO stars. The rea- 
son that these two rates do not agree more closely is not 
clear but may be related to the uncertainties concerning 
the lifetime of post-MS phases such as the RGB and HB, 
as well as dynamical considerations such as mass segre- 
gation and stellar interactions. However, even using this 
larger estimate for the turnoff rate in M15, we still find 
the formation rates implied by the thin H envelope mod- 
els to be larger than the rate at which stars are turning 
off the MS and therefore unreasonable. 

Using Figure [8] to evaluate the entire population of 
He WD candidates, 38 of these appear to have cooling 
ages < 100 Myr, with the remaining 7 either having a 
longer cooling age or being significantly redder than the 
most low-mass cooling sequence. This implies a fairly 
extreme formation rate of 380 Gyr -1 and would imply 
that between 25% to 70% of the stars turning off the 
MS are being formed into He WDs. But, many of our 
candidates are not detected in all three filters so some 
of these candidate are likely to be falsely identified, and 
we therefore can not rule out the thin H envelope models 
on this evidence alone either. However, because we find 
significant doubt as to the feasibility of the formation 
rates implied by the thin H envelope models and the 
photometric mass range for these models seems much 
less well-constrained (see previous section) we will focus 
on the thick H envelope models for the the remainder 
of this section as they appear to be the more plausible 
models. 

For thick H envelope models (Figure [9]) the 7 strongest 
candidates appear to have cooling ages down to < 1 Gyr. 
This implies a very reasonable average formation rate of 
7 Gyr -1 , and if we include the 7 bright blue gap objects 
that also seem consistent with the thick H envelope mod- 
els, we still find a reasonable average formation rate of 
14 Gyr -1 . Nine of these appear to have cooling ages < 
250 Myr implying a larger, yet not unreasonable, forma- 
tion rate of 36 Gyr -1 . In fact, this number is in excel- 
lent agreement with estimates calculated using our entire 
sample of He WD candidates in Figure [9] (see below). In 
total, 30 of the candidates plus the 7 bright blue gap ob- 
jects, appear to have cooling ages of < 1 Gyr, implying 
a formation rate of 30-37 Gyr -1 . The remaining 15 can- 



didates appear to have cooling ages of > 1.5 Gyr, which 
will result in lower average formation rate in epochs more 
than 1 Gyr ago. Although, since this cooling age is well 
below the detection limit for many of the models we likely 
have only detected a fraction of the objects with cooling 
ages between 1 Gyr and 1.5 Gyr. In summary, using the 
thick H envelope models we have obtained a formation 
rate between 7 Gyr -1 and 37 Gyr -1 depending on which 
potential He WDs are included in the sample. 

We can compare theses estimated formation rates 
with average collisional rate for RGB stars in the cen- 
tral portions of M15. W e will use the formulation in 
iBinnev fc Tremainel (|1987f ) that 

J_ = 16VW^(1 + ^). 

For our estimates we will confine the "central portion" 
to the the central 0.1' (0.3 pc). For the average stellar 
density, n, in thi s region , we w ill use the central luminos- 
ity density from lHarrisl (|1996f ) of 2.4 x 10 5 L pc -3 and 
an average M/L » 2 for the central portion of M15 from 
Ivan den Bosch et al.l (|2006| ) to get a density of 4.8 x 10 s 
stars pc -3 (assuming an average stellar mass for the cen- 
tral region of rj 1 M Q ). Using this with the velocity dis- 
persion fo r M15 a = 11 km s -1 (|Gebhardt et al.lll99i 
iDull et al"1ll997ll2003D . and the radius and mass of a red 
giant, R* = 10 R Q and M* = 0.8 M , the collision rate 
for a given red giant is approximately 1 collision per 2 
Gyr. We found about 200 RGB stars in our data that 
lie within 0.1' of the cluster center, and thus predict 100 
collisions involving an RGB star per Gyr in the central 
regions of Ml 5. 

This predicted collision rate is high enough to account 
for the upper limits on the formation rates of He WDs 
calculated above (for the thick H envelope models) under 
the assumption that He WDs are in fact formed through 
collisions involving RGB stars. However, we caution that 
our calculations here are very approximate and sensitive 
to our assumptions. 

Mass loss in a close binary system is another possi- 
ble mechanism for the formati on of He WDs. So, again 
following IStrickler et all ((20091) , we will assume that He 
WDs are formed primarily from primordial binaries and 
then estimate a lower limit on the binary fraction in the 
core of M15 by considering the formation rate of He WDs. 

Our estimated range for the He WD formation rate of 
7-37 Gyr -1 , leads to an implied binary fraction from 1% 
to 7% (using the MSTO to estimate the turnoff rate) or 
0.5% to 3% (using the HB to estimate the turnoff rate). 
These are in reasonab le agreement with the results of 
iGebhardt et alj (|1997f ) who find a binary fraction of 7% 
in M15. However, we are considering only the fraction 
of binaries that produce He WDs, which itself is only 
a fraction of the binaries in the cluster, thus our esti- 
mate should be a lower limit on the binary fraction. But, 
as shown previously in this section, collisions likely con- 
tribute to the formation of He WDs as well, and therefore 
we conclude that collisions involving RGB stars as well 
as mass transfer in close binary systems are both possible 
formation channels for He WDs in M15. 

4.4.4. Bright Blue Gap Objects as He WDs 



UV-Bright Stellar Populations in M15 



11 



Many of the bright blue gap objects are consistent with 
the thick H envelope He WD cooling models (see Fig. 
9). This, combined with the very reasonable implied He 
WD formation rates calculated with the thick H envelope 
models and the inclusion of these stars, suggests that at 
least some of the bright blue gap objects are correctly 
identified as He WDs. Many of the bright blue gap ob- 
jects are very bright and therefore appear to be quite 
young He WDs, some having cooling ages of less than 
100 Myr which may indicate that there has been an in- 
crease in the production rate of He WDs over the last 
several Myrs as a result of the dynamical evolution of 
the cluster, but is still consistent with our interpretation 
of these objects as He WDs. 

4.4.5. Cooling Ages and Implications for CO WDs 

If we use this same idea to consider the formation rate 
of CO WDs we find that 9 of the of stars we considered 
CO WD candidates have cooling ages of less than 30 Myr, 
implying a formation rate of about 300 Gyr -1 . This is a 
very small formation rate, as it is almost 2 times smaller 
than our lowest estimate for the turnoff rate for stars 
from the MS and more than 4 times smaller than the 
turnoff rate calculated using the HB. However, there are 
also 10 ambiguous WDs of which 6 appear to fall approx- 
imately in the same range for the cooling age. If these 
are added to the sample, we can calculate an implied 
formation rate of 500 Gyr^ 1 . This is a more reasonable 
formation rate, but using the turnoff rate calculated from 
the HB stars this still seems to be rather small. It is un- 
clear why this is the case, but as mentioned previously, 
there is likely some confusion in classifying CO WDs vs. 
He WDs that may account for the discrepancies in the 
formation rates. 

4.5. Nature of the Gap Objects 

In Figure Ek we have identified 60 gap objects that 
represent possible CV candidates (bright green dots and 
open circles). While this seems to be a promising set 
of candidates, when the B435 photometry is added and 
those same stars are plotted on FigureEb, most no longer 
appear to be gap objects at all but instead lie on the 
MS. Of the 60 gap objects identified in Figure EK 54 
were detected in the B435 frame as well and therefore are 
plotted on FigureEb- Two of these stars (including CV1) 
appear to be gap objects in both Figures EK & b. The 
photometry for CV1 puts it clearly in the gap region 
on both CMDs as expected. The other star mutually 
identified as a gap object, hereafter UV16 (see Table El 
and Fig. [7]), lies very near the blue edge of the MS in 
Figure Et> and almost as blue as the WD sequence in 
Figure EK- While we can make no conclusive statement 
about the nature of UV16 based on photometry alone, 
we present it as a likely CV candidate. 

Excluding CV1 and UV16, there remain 52 gap objects 
that were detected in all 3 filters. Twenty-two of these 
have NUV220 > 22 placing them in the regime where 
our NUV220 photometry is less reliable; the remaining 
30 stars with reliable photometry appear as gap objects 
in Figure Ek yet lie securely on the MS in Figure Et>- 
While this result is puzzling, our FUV140 and NUV220 
photometry is consistent with the similar results pub- 
lished by D07. In addition, these sources were inspected 



by eye to rule out spurious matches and false detections. 
Some possible explanations are discussed in the following 
sections. 

4.5.1. Cataclysmic Variable Candidates 

Since it is unclear whether the gap objects from Figure 
El are plausible CV candidates (see following sections for 
further discussion), we will focus our discussion of likely 
CV candidates to stars that appear in the gap in Fig- 
ure EH- This includes 22 objects only 6 of which were 
detected in FUV140 and therefore appear on Figure Eb- 
These 6 include CV1 and the likely CV candidate UV16 
as well 4 stars that lie either on the WD cooling sequences 
or MS in Figure Ek>. We do not consider the two that ap- 
pear on the MS as probable CVs because they show no 
UV excess in the FUV. The remaining two gap objects 
that appear on the WD cooling sequences represent much 
more probable CV candidates and have been included as 
such in Table M (UV17 & UV18) as well as shown on 
Figure[7l While both stars are dim (NUV220 > 22), they 
appear to have a UV excess in both NUV220— B435 and 
FUV140— NUV220 indicating the possible presence of an 
accretion disk. 

There are also 5 somewhat luminous gap objects (21.5 
> NUV220 > 20; 21> B435 > 19.5 - see Fig. 2a & Fig. 4a) 
that are surprisingly not detected in FUV140. We indi- 
vidually inspected these sources in the FUV140 frame and 
determined that the stars in question should be bright 
enough to be detected but are near bright stars and con- 
sequently lost in the extended PSF haloes. Hence these 
stars also remain plausible CV candidates despite the 
fact that they were not detected in our FUV140 data. 

D07 discuss several other CV candidates which they 
identified from light curves; these stars are identified in 
their paper as V7, VI 1, V15, V39, V40, and V41. Our 
photometry produced very similar results to the ones pre- 
sented in their paper although as noted previously there 
is a systematic magnitude difference in FUV140 for the 
variable stars. V15 and V40 were detected in FUV140 
but not NUV220 which is consistent with their results 
(however since they were not detected in NUV220 they 
are not included on our CMDs). We can make no further 
assertions as to the nature of these variable stars as we 
did no investigation of the variability, but we do confirm 
the photometric nature as identified in D07 as being po- 
tential CVs (see Fig. [7]for the location of V7, Vll, V39, 
and V41 on our CMDs). 

4.5.2. WD-MS Detached Binary Systems 

It is possible that at least some portion of the gap ob- 
ject population in Figure EH could be explained as WD- 
MS detached binary systems which, based on theoretical 
grounds alo ne, should be found in substantial numbers 
in GCs (see llvanova et al] 12006). We performed a very 
basic test to explore the plausibility that the photomet- 
ric properties of these gap objects could be explained as 
WD-MS detached binary systems. We fit a polynomial 
to the MS data in each CMD using a least squares fit. 
We then added the fluxes of these model MS stars to 
the fluxes from the model WD stars along the 0.6 M Q 
CO WD cooling sequence to construct a set of "observed 
fluxes for WD-MS binaries. 

In Figures [TT] & [12] we have plotted a grid where each 
line represents a specific model MS or WD star and each 
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intersection represents the resulting flux from the MS 
+ WD combination. The WD lines are solid lines la- 
beled with the appropriate temperature for the model 
WD and the MS lines are dotted lines labeled with a let- 
ter so that they can easily be followed between the two 
diagrams. For these models we only considered WDs be- 
tween 13,000 - 40,000 K and MS stars with NUV 220 mag- 
nitudes ranging from approximately 20.3 - 24.5 (models 
A- J). 

From our rudimentary model in Figures [TTJ & 1121 it is 
plausible that the observed magnitudes of the numerous 
gap objects from Figure [3k could represent a population 
WD-MS detached binaries. By examining the range of 
MS models and WD temperatures that are populated 
by gap objects in Figures [TTk and comparing them with 
Figures HTb . it can be seen that similar, though not iden- 
tical, areas of our MS- WD grid are populated in both 
CMDs. Furthermore, the majority of these objects were 
examined for variability by D07 and were not found to 
be variable. We therefore present this detached WD-MS 
binary scenario as a viable explanation for some of the 
gap population seen in Figure [3k that seems to be absent 
in Figure EJd. 

We must point out that our grid of WD-MS detached 
binary models does overlap significantly into the WD re- 
gion, specifically in the region where many of our He 
WD candidates appear (see Figure [TTk). However, we 
see that these same models do not significantly overlap 
the WD region when examining Figure [TTb . so therefore 
the addition of the optical filter allows us to more clearly 
distinguish between detached binary systems and He WD 
candidates. In addition, this overlap region in Figure [TTk 
contains a relatively large number of stars (20-25), but 
is only covered by models requiring the combination of 
WD with a very low-mass companion. So, if our He WD 
candidates arc in fact interpreted as WD-MS detached 
binaries it would indicate a strong preference for WDs to 
have low-mass companions; we find no obvious reason for 
such a preference to exist so we consider the likelihood 
that a significant amount of the He WD candidates are 
actually detached binaries to be very small. 

4.5.3. Magnetic CVs 

It is also possible that some of gap objects in Figure[5h 
could be magnetic CVs. Magnetic CVs come primarily 
in two types: polars (AM Her stars) which have no accre- 
tion disks and intermediate polars (IPs; DQ Her stars) 
which have truncated accretion disks due to most, or all, 
of the accretion in these CV systems being directed along 
magnetic field lines onto the magnetic poles of the pri- 
mary. In the field the CV population is comprised of 
abou t 25% magnetic CVs (jWickramasinghe fc Ferrarid 
2000) and it has been proposed on an observational ba- 
sis that GCs conta in a higher fraction of magnetic CVs 
than the field fe.g.. iGrindlav et aLllT9 95: Edm onds et al.l 
11999ft . In addition, llvanova et al.l (j2006ft modeled the dy- 
namical processes that lead to the formation of CVs in 
clusters and found that cluster CVs tend to have more 
massive WD primaries. Since high mass WDs have been 
show n to be more l ikely to have strong: magnetic fields 
(e.g., !Liebert|[l988t iSchmidt et al.l 12003. and references 
therein) this may explain the magnetic nature of clusters 
CVs. 

Since the accretion disk contributes strongly to the flux 



of non-magnetic CVs, polars and IPs will have different 
photometric properties than non-magnetic CVs such as 
dwarf novae (like CV1) and we argue that the missing 
accretion disk could account for the lack of NUV excess 
in Figure [3b. In FUVmo the WD dominates the flux 
while the MS companion dominates the flux in B435; the 
lack of an accretion disk leads to a significant decrease 
in the FUV140 and NUV220 flux compared to that of 
an accretion disk system and would cause the system to 
have co lors more simila r to th at of a WD-MS detached 
binary. lEdmonds et al.l (|1999l ) find that one of the CVs 
identified as a magnetic CV in their paper appears on 
the MS in the V vs. V-I CMD. Additionally, in recent 
results from GALEX observ ations of open clu ster M67, 
EU Cnc, an AM Her system (|Nair et al.ll2005ft . has been 
identified to have very similar photometric properties to 
the perplexing gap objects identified here (Hamper et al. 
2010, in preparation). 

The literature contains ideas invoking the magnetic na- 
ture of globular cluster CVs to explain observed differ- 
ences between field CVs and cluster CVs, such as the 
relatively red op tical colors and high o ptical to X-ray 
flux ratios (e.g., lEdmonds et al.l Il999t ID obrot ka et al.l 
20061 and references th erein) . lEdmonds et al.l (|2003ft and 
Dobrotka et~aH (|2006[ ) also propose that a combination 
of low accretion rates and strong magnetic moments of 
WDs in magnetic CVs may stabilize the disk such that 
outbursts would be rare, possibly explaining the low rate 
of observed outbursts in cl uster CVs (e.g., IShara et al.l 
119961 : IDobrotka et al.l 120061 and references therin). We 
therefore suggest that it is quite possible that some of 
our gap objects do represent a population of magnetic 
CVs. 

X-ray luminosity measurements would be particularly 
advantageous in constraining the nature of these objects. 
Currently 6 sets of observations of M15 from the Chan- 
dra X-ray Observatory are available in the archives. For 
a complete list of the available observations and limit- 
ing X-ray lum inosities asso c iated with each, the reader 
is referred to iHeinke et al.1 (|2009ft . Despite the gener- 
ous amount of X-ray data available for M15, the region 
of most interest for this paper, the central 20" x 20", 
is mostly dominated by two extremely bright sources 
(AC211 and X-2) that dominate most of the dimme r 
sources in this region - see Fig. 1 of IHeinke et~aH ((2009ft . 

4.5.4. Dim X-ray Source H05X-B 

Hannikai nen et al.l (|2005l ) used archival HST images 
(from the Wide Field/Planetary Camera 2) to identify 
a possible counterpart to the suspected DN H05X-B (see 
Fig. 6 in that paper). The source that we have identified 
as H05X-B is consistent with the position of the sug- 
gested counterpart for "source B" in their paper. How- 
ever, the color of this star does not appear to have any 
blue excess (see F i g. EH and Fig. [7] of this paper). 
lHannikainen et aTl ( 20051 ) identify H05X-B as a possible 
DN or qSXT, yet the color that we measure (NUV220 — 
B435 = 2.15) is redder than expected for DNe or SXTs, 
even those observed during quiesence (qSXT: Shahbaz et 
al. 2003; Edmonds et al. 2002; DNe: Bailey 1980). It is 
possible that the star that we have chosen is not the true 
counterpart. S everal other source s are w ithin the error 
circle shown in lHan nikainen ~et al.l ()2005l ) . but we exam- 
ined each of these and none seem to have particularly 
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b lue colors either. 

Hannika inen et all (|2005l) report this star as blue in 
their U-V color, so we do not rule out that the source 
identified as H05X-B could be the true counterpart. It 
is possible that the NUV220 observations in our data set 
happened to coincide with a particularly quiescent state. 
However, since the reported colors in the literature for 
quiescent DNe are bluer than those for qSXTs, it seems 
more plausible that H05X-B may be a qSXT, which is 
contra ry to the conclusion made by lHannikainen et al.l 
([20051) . 

4.6. Post-HB Stars 

We have identified 3 stars, aside from AC211, that lie 
at least 0.5 mag above the horizontal branch (see light 
blue X's in Fig. 2-4 or UV1-3 in Table 2 and Fig. 7). 
We tentatively identify these stars as being in post-HB 
evolutionary states and P-EAGB or AGB-manque can- 
didates. As these stars are only between 0.5 and 1 mag 
brighter than the HB, it is possible they have not yet 
exhausted core helium and are just beginning to evolve 
off the HB. Spectra for these stars or careful comparison 
with post-HB models is necessary to confirm the nature 
of these stars, but is beyond the scope of this paper. 
Judging solely on their relative location on the CMDs 
in Figure 21 we present UV1 as an AGB-manque candi- 
date and UV2 and UV3 as P-EAGB candidates as these 
stars app ear roughly consistent with the [Fe/H] = — 2.3 
models of iDorman et al.l (|1993[ ). 

Since P-EAGB and AGB-manque stars are consid- 
ered the progeny of the EHB, this raises the question 
of whether we should expect to see any such stars in 
M15 because we have a very small sample of EHB 
candidates. For this paper we defined the EHB as 
stars with T s ff > 20,000K, but in the low-metallicity 
models of IDorman et al.l ([1993) some stars with T c ff 
as low as 10,000K do not reach the thermal pulsation 
phase of the AGB an d therefore become P-EAGB stars. 
IDorman et al.l ([1993D find the ratio of the post-HB life- 
time to HB lifetime for AGB-manque stars to be between 
1:5 and 1:6. For P-EAGB stars this ratio varies consider- 
ably depending on luminosity, with the dimmest having 
the longest post-HB lifetime, but the ratio of lifetimes 
is usually about an order of magnitude less than that of 
AGB-manque stars (see also lBertola et al.lll995|) . 

If we treat our entire population of EHB candidates 
as true EHB stars (instead of BHks) , then we would ex- 
pect to find 1 AGB-manque star. However, it would 
seem highly unlikely to also find 2 P-EAGB stars. How- 
ever, if we include stars with T c ff > 15,000K as possible 
progenitors of P-EAGB stars, we have 16 can didate pro- 
genitors. And following IDorman et al.l (|1993| ) decreasing 
our temperature to T e g > 10,000K, we end up with 37 
candidate P-EAGB progenitors. Therefore, we do not 
rule out these stars as P-EAGB candidates, but we cau- 
tion that it is unclear whether HB stars as cool as T e g 
= 10,000K or 15,000K could produce the P-EAGB can- 
didates identified here. 

4.7. Nature of Bright Blue Gap Objects 

The nature of the objects identified as bright blue gap 
objects is very hard to determine from photometric prop- 
erties alone as they do not obviously belong to any stan- 
dard population and they populate a region of the CMD 



in which stars following canonical stellar evolution do not 
exist for any substantial period of time. From Figure 
[2}d, the 7 stars, excluding X-2, that were detected in all 
three filters and appear as bright blue gap objects on all 
CMDs (UV8 & UV10-15 in Table 2 and Fig. 7) appears 
as though they could feasibly be members of the HB or 
BS sequence; however, Figure [5^, illustrates that this is 
highly unlikely. The color of these objects makes them 
possible CV candidates, but this also seems an unlikely 
explanation because 4 of the 7 stars are brighter than 
the MSTO and several magnitudes brighter than CV1 
in all three filters. The remaining 3 have a B435 magni- 
tude similar to CV1 and the MSTO, but are significantly 
brighter in NUV220 and FUVi 40 . 

The radial distribution of the bright blue gap objects 
(Fig. [TOj) indicates that these stars are centrally concen- 
trated and therefore are likely massive binary systems. 
This, combined with the fact that X-2 (an ultra-compact 
LMXB) lies in a similar region of the CMD, leads us 
to consider that they may be close binary systems with 
some current mass transfer. If this is the case and the 
population is related to accretion disk phenomena it is 
expected that at least some of these stars should be vari- 
able. Although these stars were all included in D07, none 
were identified as variables. While this is not favorable 
for the accretion disk hypothesis, X-2 was not selected 
as variable in D07 either, as the amplitude of its vari- 
ability is less than 0.2 mag and therefore too small to be 
selected as a variable in their study. We, therefore, can 
not rule out these objects as potential mass-transfer bi- 
naries but find this explanation somewhat weak because 
none of these stars were found to show strong FUV140 
variability and none have yet been identified as sources 
of X-ray emission. 

The dimmest three bright blue gap objects (UV13, 
UV14, & UV15; Table [2]) appear to possibly be asso- 
ciated with the WD sequence, but are more than 1.5 
magnitudes brighter than the other white dwarfs in B435 
making this seem unlikely as well. Based on their posi- 
tion in the CMD alone, one could infer that the entire 
bright blue gap object population could be related to the 
BHk stars. But again, we rule this out as a likely expla- 
nation because BHk stars are not expected to be found 
more than approxim ately 1 magnitude d immer than the 
ZAHB in UV filters ([Brown et all 120011 and references 
therein); it can be seen in all the CMDs that the bright 
blue gap objects are significantly more than 1 mag. dim- 
mer than the equivalent temperature BHB stars in all 
three filters. UV8 seems to be the most likely BHk candi- 
date of the bright blue gap objects judging by its location 
in the FUV140— NUV220 CMDs, however its location in 
the NUV22o~ B435 diagrams makes its nature unclear. 

Finally, we return to the possibility that some of the 
bright blue gap objects might be young He WDs. As 
is most apparent in Figure ® these objects appear to be 
roughly consistent with the early stages of the He WD 
cooling sequence for 0.200-0.275 M He WDs in the 
thick H envelope models. In our analysis in §4.41 we in- 
cluded these stars as plausible He WD candidates when 
analyzing the cooling ages and implied formation rates 
and found that while the cooling ages implied for these 
stars require a somewhat increased production rate of He 
WDs over the last several 100 Myrs, this interpretation 
seems very reasonable. 
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We find it most likely that the majority of the bright 
blue gap objects are either young He WDs or somehow 
related to mass transfer binaries with current accretion 
disk phenomena. We find the former to be the more con- 
vincing explanation but do not rule out that the bright 
gap population may encompass two or more physically 
distinct populations, especially since it does include X-2. 

5. DISCUSSION 

5.1. Gap Population 

The filter combination of FUV140 an d NUV220 was 
very successful in identifying members of the gap zone, 
which are potential CVs as pointed out by D07. How- 
ever, the addition of the B435 filter has made the nature 
of these stars less clear. Because the majority of these 
stars appear on the MS in NUV220 — B435 CMDs, it seems 
unlikely that they are standard DNe CVs. We find the 
most likely explanation to be that this group is composed 
of a combination of magnetic CVs and WD-MS detached 
binaries. There is no analog in the literature for such 
a significant population of FUV-bright objects that ap- 
pear MS-like in a color such as NUV220— B435; so we can 
draw only weak conclusions about the gap population we 
have identified in Figure [3^ without the addition of Ha 
imaging to seek out emission line sources or spectroscopic 
follow up. 

5.2. M15 as a Blue Hook Cluster 

We have identified 5 of the EHB stars as potential 
BHk stars. If we are correct in our classification of these 
stars as BHk stars, it may provide an important data 
point in understand ing the properties of BHk clusters. 
iDieball eT al. (2009) investigated trends among clusters 
containing BHk stars with their strongest result being 
that BHk stars seem to exist exclusively in the most mas- 
sive clusters. They note that this might be attributed 
to a natural bias toward massive clusters, since intrinsi- 
cally rare stars are necessarily more likely to be found in 
larger samples. They also found weaker correlations be- 
tween the presence of BHk stars and other cluster prop- 
erties such as concentration parameter, core radius, and 
relaxation time, but considered the most significant cor- 
relation to be that with cluster mass. Ml 5 was included 
in their sample, but it was considered a non-BHk clus- 
ter as it was only known to have one BHk star in the 
outer region (see §3.4.3|) . Had Ml 5 been included as a 
BHk cluster in that study it would have had a strong 
lever arm on the results, because, of the clusters in their 
study that had 4 or more BHk stars, M15 has the low- 
est mass, highest central density, shortest core relaxation 
time, smallest core radius, and highest concentration pa- 
rameter. However, it is unclear whether M15 would be 
considered a BHk cluster even with all the candidates 
we have ide ntified here in c luded since all the "BHk clus- 
ters" in the IDieball et all (j2009D sample had 10 or more 
BHk stars. Yet, we have only searched a small region in 
M15, so there may be more BHk stars that we have not 
identified. Nevertheless, irrespective of how we define a 
BHk cluster, it is clear that a BHk population in M15 
provides an important constraint in understanding the 
origin of BHk stars. 

Since we feel that the existence of a BHk population 
most likely weighs in favor of the late He-flash scenario, 



as discussed in £14.31 it should be noted that within this 
scenari o the expected mass ra n ge for BHk stars is quite 
small. iMiller Bertolami et al.l ((2008) calculate that the 
expected mass of a low-metallicity remnant capable of 
becoming a BHk star via a late He-flash lies in a very 
narrow range between 0.48 - 0.50 M Q (see that paper for 
more details) . It is expected that there would be a spread 
in the total mass loss from the progenitors of these BHk 
stars thus resulting in a comparable number of stars with 
post-MS ages similar to BHk stars (< 100 Myr) that just 
"miss" becoming BHk stars and end up as either massive 
He WDs (which just avoid He-core ignition) or low-mass 
CO WDs (which arrive as such following a phase as an 
EHB or BHk star). It is unclear whether M15 contains 
such a population. There are a significant number of 
potential WDs that may be either low-mass CO WDs or 
massive He WDs in the appropriate age range (see Fig. 
[8] & EJ, but many of these candidates have ended up 
classified as "Ambiguous WDs" (Tabic [3]) or were only 
detected in two frames due to the intrinsic photometric 
uncertainty associated with such dim stars. So, without 
a more precise determination of both the mass and age 
of our WD candidates, we can make no further claim to 
the existence of such a population other than to note its 
expected presence. 

6. SUMMARY & CONCLUSIONS 

We have presented a photometric identification and 
analysis for several UV-bright populations in the cen- 
tral region of the post-core collapse globular cluster M15. 
We additionally have included photometry for 4 previ- 
ously identified X-ray sources: M15 CV1, AC211, M15 
X- 2, and H05X-B. Our work has elaborated on the work 
of IDieball et al.l (|2007f l. who analyzed the FUVi 40 and 
NUV220 images. We reanalyzed these images and added 
the B435 filter and the NUV220 - B435 color which has 
added further insight into the nature of the populations 
discussed here. The UV-bright populations we have ana- 
lyzed include many stages of non-canonical stellar evolu- 
tion including blue stragglers, extreme horizontal branch 
stars, blue hook stars, helium-core white dwarfs, and cat- 
aclysmic variable candidates. 

We have selected 53 blue straggler candidates which 
display a clear central concentration as expected for 
BSs. Since the expected zone of avoidance is beyond 
the field of view for our images we are unable to inves- 
tigate whether the BS population displays the double 
peaked bimodal distribution that has been discussed in 
the literature for several other clusters. We found 60 
CV candidates populating the gap between the MS and 
WD region, consistent with the findings of D07, however 
upon inclusion of the B435 filter we found that many of 
these stars do not display expected CV colors, but in- 
stead appear as MS stars. Thus, we suggest that these 
gap objects may be magnetic CVs with truncated or ab- 
sent accretion disks or possibly detached WD-MS bina- 
ries. However, we also find three gap objects that we 
consider very likely CV candidates (UV16-18) as they 
display colors similar to what is expected for a CV. 

We have used a ZAHB model to select 6 extreme hor- 
izontal branch candidates, 5 of which appear to be sub- 
luminous in the UV and therefore better candidates for 
blue hook stars. Due to the existence of these candi- 
dates, in addition to the one previously identified BHk 
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star, we suggest that M15 be considered a BHk cluster 
for future studies on clusters containing BHk stars as it 
may provide important constraints on how these stars 
are formed. We also identify three stars that represent 
plausible post-HB stars that may be AGB-manque or 
post-early AGB stars, the progeny of EHB stars. Our 
identification of these post-HB stars however is very pre- 
liminary as their photometric properties seem consistent 
with AGB-manque and P-EAGB states, yet statistically 
it seems improbable to have detected three stars in these 
relatively short lived stages of evolution. 

Additionally, we uncovered a population of "bright 
blue gap objects" for which there is no obvious ana- 
log in the literature. We consider these stars to most 
likely be young He WDs, but they could be related to 
the BHk population or accretion disk binary systems. 
UV10-UV15 seem to be the most plausible candidates 
for He WDs. There is no obvious population to which 
UV8 belongs but it seems to be related to the HB and 
is a possible BHk candidate. In addition to the bright 
blue gap objects, we have identified a significant pop- 



ulation of candidate He WDs which we analyzed using 
model He WD cooling sequences. We have probed quite 
deep into the He WD sequence and find 7 strong He WD 
candidates detected in all three images; an additional 38 
candidates were detected in only two images. This is the 
first strong evidence of the existence of a He WD popula- 
tion in M15. We analyzed both thin and thick H envelope 
models and based on the cooling ages and overall fit to 
the population we find the thick H envelope models to be 
better fit to our data. The formation rates suggested by 
the thin H envelope models are unreasonable unless we 
have misidentified several of our candidate He WDs. The 
formation rates implied by the thick H envelope models 
range from 7-37 He WDs produced each Gyr. The col- 
lision timescale of red giants in the core of M15, is high 
enough that it is possible that collisions account for a 
significant fraction of the He WDs. Furthermore, the im- 
plied lower limit on the binary fraction calculated from 
these formation rates suggest that close binary systems 
are also likely to contribute to the formation of He WDs 
in M15. 
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TABLE 1 

Previously Identified X-ray Binary Systems 



Star 


R.A. 


Dec. 


B435 


NUV220 


FUV140 


AC211 a 


21:29:58.323 


+12:10:01.94 


15.23 


14.49 


13.91 


X-2 6 


21:29:58.142 


+12:10:01.52 


19.51 


17.98 


16.96 


CV1 C 


21:29:58.357 


+12:10:00.33 


20.42 


20.28 


21.06 


HX05-B d 


21:29:58.323 


+12:10:11.69 


21.19 


23.33 





Ref erences. — lAuriere et"aLl (|1984f ): b IDieball et all (120051 ); IShara et all (120041 ); IHannikainen et al.1 
(2005) 



TABLE 2 
UV-Bright Stars of Interest 



Star ID 


XF220W 


YF220W 


B435 


NUV220 


FUV 140 


Candidate Classification 


UV1 


233.98 


93.13 


16.14 


15.57 


14.85 


AGB-manquc 


UV2 


341.12 


920.76 


14.61 


16.22 


16.82 


P-EAGB 


UV3 


579.31 


500.38 


14.74 


16.42 


17.74 


P-EAGB 


UV4 


773.86 


722.57 


17.88 


16.70 


15.67 


BHk 


UV5 


65.31 


1016.84 


18.07 


16.89 


15.85 


BHk 


UV6 


672.66 


595.49 


18.40 


17.04 


15.97 


BHk 


UV7 


932.36 


976.88 


18.40 


17.20 


16.16 


BHk 


UV8 


688.90 


696.32 


17.74 


17.20 


16.47 


Unknown 


UV9 


1028.58 


1072.23 




17.42 


16.08 


BHk 


UV10 


636.65 


635.43 


18.70 


18.31 


17.45 


He WD 


UV11 


548.94 


603.81 


20.10 


18.43 


17.16 


He WD 


UV12 


630.50 


602.71 


18.58 


18.55 


17.83 


He WD 


UV13 


629.48 


698.07 


18.82 


18.69 


17.86 


He WD 


UV14 


582.81 


643.68 


20.47 


19.01 


17.85 


He WD 


UV15 


412.42 


687.64 


20.07 


19.01 


18.06 


He WD 


UV16 


209.254 


608.48 


20.24 


20.95 


20.12 


CV 


UV17 


626.60* 


369.06* 


22.33 


23.24 


21.81 


CV 


UV18 


212.69* 


980.50* 


23.42 


24.25 


23.68 


CV 



Note. — The coordinates listed as XF220W & YF220W are x, y coordinates for the sources in the archival NUV220 image 
"j8si08010_drz.fits. " The coordinates for UV17 & UV18 are marked with asterisks as these sources are both very dim and 
therefore difficult to identify in the this frame, but are clearly present in the fully combined images we used as our master 
images. 



TABLE 3 

Summary of WD candidates 





Total Candidates 


Detected in All Three Frames 


All WDs 1 


73 


20 


CO WDs 


11 


5 


He WDs 


45 


7 


Ambiguous WDs 2 


10 


5 


Possible WDs 3 


5 


3 


D07 Variables 


2 




Bright He WDs 4 


7 


7 



Note. — 1 "All WDs" does not include the bright blue gap objects, as this population is dis- 
connected from the canonical WD region. 2 "Amibguous WDs" are strong WD candidates which lie 
near the transition between the CO and He WD cooling sequences, making their identification as a 
CO vs. He WD very difficult. 3 "Possible WDs" are candidates with significant photometric error or 
inconsistencies such that their nature is unclear. 4 "Bright He WDs" are bright blue gap objects that 
appear consistent with He WD cooling sequences. 
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Fig. 2. — Comparison of NUV220 — B435 and FUV140— NUV220 CMDs. Stellar populations are distinguished and color-coded according 
to their positions in 2a, the NUV220 - B435 CMD. See key upper righthand corner of 2b for symbol explanation. Theoretical WD cooling 
sequences for CO WDs (solid lines) ranging in mass from 0.45 - 1.10 Mq and thick H envelope He WD s (d ashe d lines) ranging in mass from 
0.175 - 0.45 Mq, as well as two ZAHBs (solid dark red line and dark blue lines) are included; see £|4,ll and l3.4l for details. All variable stars 
were identified as such by D07. The "dim gap objects" (NUV > 22) have been plotted as a separate group so they can be distinguished 
on all CMDS, the reason for this differentiation is described in i|4.2.3l Squares mark optical counterparts of four X-ray sources discussed 
in the text. 




Fig. 3. — Same as Figure [2] except populations distinguished and color-coded according to their positions on the FUV140— NUV220 (3a) 
CMD. See key in upper righthand corner of 3a for symbol explanation. The RGB has not been distinguished from the MS in these figures 
as it is very difficult to the distinguish the RGB in 3a, as discussed in the text. 
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(») 




main sequence 

* red giant branch 

■ asympt. giant branch 
horizontal branch 

■ post horiz. branch 

n extreme horiz. branch 
^possible blue hook 

* bright blue gap objects 
T blue stragglers 

A white dwarfs _ 
gap objects 
dim gap objects 
A variable stars 



(b) 





• not detected in B 



Fig. 4. — Similar to Figures [2] & [3] but with the magnitude axis being NUV220 for both CMDs. Populations distinguished and color-coded 
according to their position in Figure [2jt. See key in upper righthand corner of 4a for symbol explanations. 




Fig. 5. — The location of the suspected counterpart to H05X-B in B435, NUV22O) and FUV140 respectively. 






Fig. 7. — CMD showing the stars from Tables [T] & [2] (red squares) as well as CV candidate variable stars from D07 (purple triangles). 
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Fig. 8. — WD candidates and model cooling sequences for CO WDs (blue solid lines) and thin H envelope He WDs (purple dot-dash 
lines). Masses (in M ) for models, from left to right - CO WDs: 1.10, 1.00, 0.90, 0.80, 0.70, 0.60, 0.50, 0.45; thin H envelope He WDs: 
0.45, 0.35, 0.30, 0.25, 0.20, 0.175. Cooling ages are marked along the cooling curves and indicated in the key located in the upper right. The 
cooling curve for the fiducial 0.6 Mq CO WD has been plotted as a thicker line for orientation purposes. Filled circles are WD candidates 
that were detected in all three frames; open circles were only det ected in two frames; larger red filled circles were detected in all three 
frames and represent the 7 strongest He WD candidates (see ij4.4t ; and filled triangles are stars identified as varia bles in D0 7. The grey 
asterisks are bright blue gap objects that we consider possible He WD candidates based on the curves in Figure l9l f i|4. ll fc 1 1.71 1. Error bars 
shown were calculated by the program ALLSTAR in DAOPHOT II (see Stetson et al. 1990). The error bars in the upper portion of the 
figure show the lo~ error for the cooling curves due to the uncertainty in distance and reddening. 
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Fig. 9. — WD candidates and model cooling sequences for CO WDs (blue solid lines) and thick H envelope He WDs (green dashed lines). 
Symbols for the WD candidates, bright blue gap objects, and CO WD curves are the same as in Fig. [8] Masses (in Mq) for models from 
left to right - CO WDs: 1.10, 1.00, 0.90, 0.80, 0.70, 0.60, 0.50, 0.45; thick H envelope He WDs: 0.45, 0.40, 0.35, 0.30, 0.275, 0.25, 0.225, 
0.20, 0.175. Cooling ages for thick envelope He WD curves are indicated in the key located in the upper right. 
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Distance From Cluster Center (arcsec) 



Fig. 10. — Cumulative radial distributions for selected stellar populations. The populations used for the left panel were selected from Fig. 
2a and the populations used for the right panel were selected from Fig. 3a (i.e. "BS Candidates" on left panel represents the distribution 
of those plotted as inverted blue triangles on Fig. 2 and "BS candidates" on the right panel represents the distribution of those plotted as 
inverted blue triangles in Fig. 3). See key in lower righthand corner of right panel. The distribution for "Gap Objects" includes all objects 
identified as gap objects (no magnitude limit). 
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Fig. 11. — Gap objects from Figure^ (green dots) are plotted with grids from our MS- WD detached binary system models. Gap objects 
that were not detected in B435 are plotted as green inverted triangles and stars classified as WDs are plotted a s grey pinched triangles. 
Intersections are labeled by the effective temperature of the WD and a letter representing the MS model (see £|4,5.2jl ; each intersection 
represents the resultant flux for the combination of the MS and a 0.6 Mq CO WD star. MS models range from NUV220 ~ 20.3 - 24.5 
(A- J). 
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Fig. 12. — Color-color diagram for stars detected in all three filters; symbols as in Fig. 1111 For reference, MS stars appear as a clump on 
the right hand side, HB stars appear at the top above the model WD-MS binaries, and the BSs form a sequence near the line corresponding 
to WD-MS binaries that include the MS model labeled F. 



